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Abstract. There are two types of particles interacting on a homoge- 
neous tree of degree d + 1. The particles of the first type colonize the 
empty space with exponential rate 1, but cannot take over the vertices 
that are occupied by the second type. The particles of the second type 
spread with exponential rate A. They colonize the neighboring vertices 
that are either vacant or occupied by the representatives of the oppo- 
site type, and annihilate the particles of the type 1 as they reach them. 
There exists a critical value A c = (2d — 1) + \J (2d — l) 2 — 1 such that 
the first type survives with positive probability for A < A c , and dies out 
with probability one for A > A c . We also find the growth profile which 
characterizes the rate of growth of the type 1 in the space-time on the 
event of survival. 



1. Introduction 

We consider a model of a predator-prey type which we call the Escape 
model. There are two entities growing on the vertices of a homogeneous tree 
Td = T of degree d+1. The entities may be thought of as biological species, 
political parties or manufacturers competing on a market. The second entity 
dominates the first in the sense that the representatives of the second entity 
can take over the vertices occupied by the representatives of the first entity 
but not vice versa. It is also assumed that the second entity grows faster. We 
are interested in the possibility of the long-term coexistence of the species 
which occurs when the first species survives. 

At each time t each site of the tree is occupied by at most one repre- 
sentative of either the two entities. We refer to the representatives of the 
entities as particles of types 1 and 2 respectively. If a site is not occupied, it 
is said to be vacant. The dynamics of the process is specified by exponential 
transition rates. A vacant site gets colonized by a particle of type 1 with 
exponential rate equal to the number of nearest neighbors of type 1. The 
sites that are either vacant or occupied by type 1 flip to 2 with the rate 
A > 1 times the number of neighbors of type 2. If a vertex is occupied by a 
particle of type 2, the particle stays at the vertex forever. We assume that 
at time zero there are finitely many particles of each type. 

l 
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The sets of sites occupied by particles of types 1 and 2 at time t > 
are denoted by A(t) and B(t) respectively. The event r\t>o{A(t) ^ 0} is 
referred to as the event of survival of type 1. If at time zero the cluster ^4(0) 
is surrounded by B(0), then all particles of type 1 eventually die out with 
probability one. Thus, we are only interested in the initial configurations 
for which there exist a vertex x and an infinite geodesic segment 7 Xj0 o such 
that x € A(0) and "f x ,oc H -B(O) = 0. All such configurations are referred to 
as non-trivial configurations. The question of interest is for what values of 
A type 1 survives with positive probability in the long term. 

Theorem 1. There exists a critical value A c = (2d — 1) + \J (2d — l) 2 — 1 
such that, for all A G (1, A c ) and for all non-trivial finite configurations 
(A(0),B(0)), type 1 survives with positive probability. For all A € (A c ,oo), 
type 1 dies out with probability one. 

For c > denote by M n (n/c) the number of vertices x S A(n/c) at 
distance n from the root p. From the results of Sections 0] and 03 it follows 
that 

lim-log(EM n (n/c)) = -g(c). 
n 

where 

( (A/c-log(A/c) - 1) -logd : 0<c<l 

g(c) = l (A/c-log(A/c)-l) + (l/c-log(l/e)-l)-logd : Kc<A 
[ (1/c - log(l/c) - 1) -logd : c>A 

The function g(c) is referred to as the growth profile of type 1. The growth 
profile was introduced in Lalley |3I (in a slightly different form) to study the 
weakly supercritical contact process on a homogeneous tree. The function 
g(c) has a unique minimum at cq, is strictly decreasing on (0, Co) and strictly 
increasing on (co, oo). For all A G (1, A c ), we have g(co) < 0, and let r\ and 
r2, with < n < r2, be the two solutions of g(c) = 0. 

Theorem 2. Lei A € (1, A c ). .For every e > and a// /arye i/ie particles of 
type 1 are concentrated in the annulus of radii (r\ — e)t and (r2 + e)t centered 
at the root. For every c G (^1,^2); almost surely on the event of survival of 
type 1, 

lim- log (M n (n/c)) = -g(c). 
n 

To investigate the Escape model, we consider the simplest growth model, 
the Richardson model, presented in Sectional 

2. Preliminaries. 

2.1. A homogeneous tree. A homogeneous tree 7^ = T of degree d+ 1 is 
an infinite tree such that every vertex has exactly d + 1 nearest neighbors. A 
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distinguished vertex is called root and denoted by p. For every two vertices 
x and y of the tree, denote by \x, y\ the number of edges in the shortest path 
from x to y (the path without loops). If x is the root, then we simply write 
\y\. Note that |-, -| is a metric on T. 

For every vertex x G T, a geodesic segment ^y x ,oo is an infinite path in T 
beginning at x and having no loops. Define by T + (x), the set of all vertices 
y such that the shortest path connecting y with the root p goes through x. 
Consider a geodesic segment 7^00 such that 7^00 G 1+(x). Given an integer 
m, consider also a sequence of vertices {yk)k>o on the geodesic segment j x ,oo 
such that \yk\ = mk + \x\ for all integers k > 0. For every k > 1, y^-i is 
called the m-predecessor of in T. 

Let D(x,r) = {y G T : |x,y| < r} be the closed disk of radius r G (0, 00) 
centered at x, and let C(x,r) = {y G T : \x,y\ = r} be the circumference 
of that disk. If x is the root, then we write D r and C r . Note that for all 
integers n > 1, the number of vertices in C n is (d + l)d n . 

2.2. Construction of the Escape process. The Escape process can be 
built using a percolation structure as follows. For each ordered pair of neigh- 
boring vertices x and y in T, define two independent Poisson processes with 
respective rates 1 and A — 1, and respective occurrence times {Tn' y : n > 1} 
and {Un' y ■ n > 1}. Set Tq ,y = and C/q' 27 = 0, and make these Poisson pro- 
cesses independent from pair to pair. Consider T x R + . Arrows are drawn 
from x to y at the occurrence times Tn' v and Un V ■ We say that there is a 
directed path in T x M + from (xo, so) to s n +i) if there is a sequence of 
times sq < si < .. < s n+ i and sequence of vertices xo,x\, ..,x n so that for 
each j, 1 < j < n, there is an arrow from Xj-i to Xj at time Sj. A type 1 
pai/i is a directed path that uses only arrows generated by the Poisson pro- 
cesses T n . Fix an initial configuration A(0),B(0) and erase all arrows that 
lie only on paths that begin at points (x, 0) such that x g" (j4(0)UB(0)). For 
every vertex y we say that y G B(t) if and only if there is a vertex z 6 B(0) 
and directed path from (z, 0) to (y, t) in the modified percolation structure. 
Define A(t) to be the set of vertices y such that y £ B{t) and there is a type 
1 path (in the modified percolation structure) that ends at (y, t) and starts 
at (z, 0) for some z G A(0). 



3. The Richardson model on a homogeneous tree 

The Richardson process on T with parameter A > is a continuous time 
Markov process R(t) on the set of finite subsets of T. We say that a vertex x 
is infected (or occupied) at time t if x G R(t), and is vacant otherwise. The 
process develops according to the following rules: if a vertex gets infected, 
it stays infected forever and starts infecting unoccupied neighboring vertices 
with rate A, i.e. that the infection times have exponential distributions 
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with parameter A. The infection times are all independent. Consequently, 
a vacant site becomes infected with the rate 

A • (number of infected neighbors). 

The model with parameter A = c can be obtained from the model with 
parameter A = 1 by running the later process at speed c. Therefore we treat 
just the case A = 1. 

Without loss of generality consider the initial configuration where at time 
zero the only occupied site is the root R(0) = {p}. The main questions were 
how fast the infected set grows and what limiting shape the infected set has. 
For the tree T the number of vertices in the disk D n grows exponentially with 
n (whereas for Z d it has polynomial growth). Consequently, the behavior 
of the Richardson process on T is different from the behavior of the model 
on the integer lattice (see 0j and U). The infected region R(t) still grows 
linearly with time, but there are constants a and 6, with a < b, such that 
as t goes to infinity we can classify two subregions: a completely infected 
subregion, having approximately the shape of a ball Dt a , and a partially 
infected subregion, having approximately the shape of a ring \ Dt a - 
Consequently, there are two speeds: the speed of invasion b, indicating how 
fast the infection spreads, and the speed of occupation a, governing the rate 
of growth of the region that is completely covered by the infection. 

Proposition 1. Let d > 2 be an integer. Let 

(1) /(c) = --log--l-logd 

c c 

for c € (0, oo), and let a and b, with < a < 1 < b, be the two roots of the 
equation 

f(c) = 0. 

Then a = sup{a'} and b = inf{6'}, where the sup and inf are taken over all 
a' and b' satisfying 

P[3 a random r < oo such that, Vi>r, Dt a > C R(t) C Dtv\ = 1. 

As d^oo, we have a(d) log d — > 1 and b(d)/d — > e. 



Let N n (t) be the number of vertices at distance n from the root that are 
infected at time t, and let F n (t) be the number of vertices in C n that are 
not infected at time t. We compute asymptotic values of N n (t) and F n (t) 
as n goes to infinity and t = n/c for different values of c. 

Proposition 2. For all c G (1,6), 

(2) lim - log (N n (n/c)) = -/(c) > a.s.. 

n 

For all c € (a, 1), 

(3) lim - log (F n (n/c)) = -/(c) > a.s.. 
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(4) P[3r < oo such that Vt > r, A(o-e) C -R(t) C A(6+e)] = 1j 



(5) 



P[3t < oo such that W > r, C A(b-e)] = 0, 



(6) 



P[3r < oo such that Vi > r, A(a+e) c = °- 



Proof of equation (fljl. For every vertex x on the tree, let T(x) be the (ran- 
dom) time at which the vertex gets infected. Consider an arbitrary vertex 
at distance n from the root, and denote it by x n . Fix c > 1. Then we have 
the following estimate: 



n 



P [x n € R(n/c)} = P T{x n ) < - 
' L c 

where 



exp 



-n ( log 1 ) + o(n^ 



o(n)/n — > as n — > oo. 

The estimate follows from Cramer's theorem for i.i.d. random variables 
(see P]) and from the fact that T(x n ) is distributed as a sum of n i.i.d. 
exponentials with parameter 1. Observe that 

EN n (n/c) = (the total number of vertices in C n ) ■ P(x n E R(n/c)) 



d+1 
d 



d n exp 



1,1 

-n | log 1 I + o(n) 

c c 



(7) 



exp{— nf(c) + o(n)} 



where / was defined in lfl|). 
Similarly, for c < 1 we have: 



P[x n ^ i?(n/c)] = P T(x n ) > - 

c 



n 



exp 



-n ( log 1 ) + o(n) 



Hence 

(8) E F n (n/c) = exp{-n/(c) + o(n)}. 

Observe that / is strictly decreasing on (0, 1) and strictly increasing on 
(1, oo), with unique minimum at c = 1. Moreover /(l) < and /(0+) = 
/(oo— ) = oo. Thus, there are just two roots a and b of the equation /(c) = 0, 
such that a < 1 < b. Hence, for every e > 0, we have that 



3x G C n : T(x) < 



n 
b + l 



<EiV n (n/(6 + e)). 
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Since f(b + e) > 0, by J2f the upper bound decays exponentially with n. 
Consequently the left side is summable, and by Borel-Cantelli lemma we 
have 



3N < oo : Vn > N and Vx G C n , Tlx) > 



n 



b + e 



1. 



Finally, observe that this is equivalent to 

(9) P [3r < oo : Vt > r, R(t) C £> t(He) ] = 1 

(the events are identical). 

Analogously, for all e G (0, a), 

(10) P [3t < oo : Vt > r, D t (a-e) C = 1. 
To prove note that 

n 



3x G C n : T(x) > 



<EF n (n/(a-e)), 



and it decays exponentially by Applying Borel-Cantelli lemma one more 
time, we get 



< oo : Vn > N and Vx G C n , T(x) < 



1 



which implies (|1U|) . Obviously, © and (|1U|) together are equivalent to @. 

Equations © and © are direct consequences of Proposition |21 As a 
heuristic argument, note that © implies that, for every c G (1, b), EN n (n/c) 
grows exponentially, and, similarly, © implies that for every c G (a, 1), ~EF n (n/c) 
grows exponentially. These observations suggest that © and © should be 
true. 

The limits 

lim a{d) log d = 1, 
lim 6(d) /d = e 

immediately follow from the fact that a and 6 are the roots of 

1 - log d = 0. 



1 i 1 

- - log - 

c c 



The phenomenon is easily anticipated. Since for larger d's there are more 
directions for the infection to spread around, it is natural that the invasion 
speed is strictly increasing with d. For the same reason, the occupation 
speed decreases to zero (the number of vertices in D n grows unboundedly 
with d). □ 



Proof of Proposition We only prove © . The proof of © is identical. 
First we claim that for any e > 

lim sup - log (N n (n/c)) < -/(c) + e 
n 
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almost surely. By Markov's inequality and (jJJ), 

P [N n (n/c) > exp {n(-/(c) + e)}] < 



-log(N n (n/c)) > -/(c) + e 
n 



< E [JV n (n/c)] • exp {nf(c) — ne} = exp {— ne + o(n)} , 

and the claim follows by Borel-Cantelli lemma. 

Let a = -. To finish the proof, it is enough to show that, for an arbitrarily 
small e > 0, with probability 1 

(11) liminf - log (NJan)) > -f (1/a) - e 

n 

Recall that 1 < c < b, so / (1/a) < 0. Assume that e is small enough to 
satisfy / (1/a) + e < 0. To make the notation less complicated, let 

H = n( e ) = exp {- (/ (1/a) + e)} . 

By the continuity of /, there exists an e% > such that 

/(l/(a-2ei))</(l/a)+e<0. 

Then, for every fixed w > 1, there exists an integer m > large enough such 
that 

EN m ((a - 2ei)m) > w ■ exp {-m (/ (1/a) + e)} 

(12) = w ■ fi m > 1. 

Choose m to satisfy ()12l) . Fix an arbitrary vertex x of the tree and consider 
a geodesic segment 7^00 £ T+(x). Consider a sequence of vertices (yi)i>o 
on the geodesic segment 7^,00 such that \yi\ = mi + \x\ for all integers 
i > 0. Note that infection times T{yi) are increasing in i. For every pair 
of non-negative integers m and 722 such that m < rt2, say that y n2 is an 
m-descendant of y ni if, for all integers i G [m, 71-2), 

- T{yi) < (a - 2ei)m. 

Define 

Zk{x) = {z £ Ck m+ \ x \ : z is an m-descendant of x} , 
Zk(x) = cardinality oiZ^{x). 

Note that (Zfc(x))fc>o is a Galton- Watson process with mean offspring num- 
ber 

(13) EZi(x) > w-n m > 1. 

Claim 1. For every vertex x, almost surely on the event of (Zfc(x)) fc>0 
survival, there exists a (random) K < +00 such that, for all k > K, 

(14) AWb|((" " ei)(mk + |x|)) > Z k (x) > fi mk +^ ■ d m . 
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Proof. The first inequality follows from the fact that for sufficiently large 
fc's 

(a — ei)(mk + \x\) > (a — 2e\)mk + T{x). 

To obtain the second inequality, observe that, for all large enough k, w k > 
d m jj} x \, and hence, EZ k (x) > (/x m • w) k > ^ mk +\ x \ . d m . A standard theo- 
rem from the elementary theory of Galton- Watson processes states that if 
EZi(x) > 1 and the variance of Z\{x) is finite, then on the event of survival 

Zj? /(E z[ x ^) k converges almost surely to a positive random variable. Thus, 
the second inequality is obtained by direct application of (|T3|) . □ 

Therefore, on the event of non-extinction of (Zk(x)) k>0 , (jTI|) is true and 
implies for n's of the form mk + \x\. To establish the result for all 
positive integers, consider mk + \x\ < n < m(k + 1) + \x\. Observe that K 
might be also chosen large enough that, for all k > K, we have 

(a — ei)(m(k + 1) + \x\ ) < a(mk + \x\ ). 

By ()14|) at time (a — ei)(m(k + 1) + |x|) we have at least 

^m(k+l)+\x\ _ ^rra 

infected vertices on level m(k + 1) + \x\. Since each particle can generate 
at most d offspring, it follows that, for each integer n satisfying mk + \x\ < 
n < m(k + 1) + \x\, there are at least 

^m(k+l)+\x\ ( jn—(mk+\x\) 

infected predecessors in C n at time (a — e\){m(k + 1) + Obviously, 
N n (an) > N n (a(mk + \x\)) 

> N n ((a - ei)(m(A; + 1) + \x\)) > fl m ( k + 1 )+\ x \ > ^. 

Therefore we proved that, almost surely on the event of the survival of 
(Zfc(x))fc>o, is true. To show that holds with probability one, 
observe that, for each integer j > 0, there are (d + X)d m3 ~ Galton- Watson 
processes (Zk(x))k>o with \x\ = mj. Let Sj be the event of non-extinction 
for at least one of the processes. Since the Galton- Watson processes are 
independent, the probability of Sj tends to 1 as j tends to +oo. Moreover, 
for each j > 0, Sj C Sj+\ which guarantees the almost sure result. □ 

4. TWO INDEPENDENT RICHARDSON MODELS. 

To motivate the proof of Theorem^consider two independent Richardson 
processes Ri(t) and R\(t) with respective rates 1 and A > 1. The processes 
have initial configurations -Ri(O) = R\(0) = {p}, and are built on a homo- 
geneous tree T. 

For c £ (0, oo), estimate the expected number of vertices in (R\(n/c) \ 
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R\{n/c)) n C n , that is the number of vertices at distance n from the root 
that are occupied by R\ but not by R\ at time n/c. Let x n be a vertex 
with \x n \ = n, u n (n/c) = P [x n G {Ri(n/c) \ R\(n/c)}], and V n (n/c) = 
# {x n ■ x n G R\{n/c) \ R\(n/c)}. 
Case 1: For every c G (0, 1], 

u n (n/c) = P [x„ G" i?A(«/c)] P [s„ G J?i(n/c)] = 

= exp < — n log 1 + o(n, 

{ \c c J 

Recall that the number of vertices in C n is ^-d n and define 
9i(c) = - log^ - 1^ - logd. 

Then, for all c G (0,1], 

EV n (n/c) = exp{-n • 51(c) + o(n)} . 

Case 2: For every c G (1, A) 

u n (n/c) = P[x n <£ R\{n/c)\ P [x n G R\{n/c)} = 

exp I — n ( log 1 J + o(n) 1 exp { — n ( log 1 ) + o(n) 



c c j J L\ c c 

Let 

52(c) = - log^ - 1^ + Q - log^ - 1^ - logd. 

Then, for all c G (1, A), 

EV n (n/c) = exp{-n • 32(c) + o(n)} . 
Case 3: For every c G [A, 00), 

u n (n/c) = exp < — n ( log 1 ] + o(n) 

{ \c c ) 

Thus, 

EV n (n/c) = exp{-n • 33(c) + o(n)} 

where 33(c) = Q - log \ - l) - logd = /(c). 

Define function 3(c) on (0, 00) by combining 31(c), 32(c) and 33(c) on their 
domains. Note that for every c > and non-negative integers m and n 

u m {m/c)u n {n/c) < u m+n ((m + n)/c). 

Hence 

(15) EV n {n/c) < exp {-n ■ 3(c)}. 

The function 3(c) is continuosly differentiable, strictly decreasing on 

(0, (A + l)/2) and strictly increasing on ((A + l)/2, 00) with the unique 

minimum at co(A) = A±l. Furthermore, 3a (co) = log . It follows that 

3(c ) > if A > (2d - 1) + v / (2d- l) 2 - 1, 
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g(co) < if 1< A < (2d - 1) + ^(2d - l) 2 - 1. 



Proposition 3. Fix A > A c = (2d — 1) + (2d — l) 2 — 1. Almost surely, 
for all sufficiently large t, 

R x (t) c R x (t). 

Proof. According to Proposition ^ for any e\ > and all large t, 

Thus, we are only interested to see what happens in the region ^(oA-ei)' 
Note that by (|15[) . for every e > and large n, 

E [# {x n : x„ G Rx(n/c) \ R\(n/c - 1)}] < exp {-n • 5(c) + ne} . 

For large t > and integers i > 1, let q = i/t. By Markov's inequality, 

00 

(i? x (i) \ J2 A (t - 1)) n D c t{aX _ ei) + 0] < ^ exp {-i • g{a) + ie} < 

i=[t(aA— ei)] 

00 

< ^ exp{-i5f(c ) +ie} < Cexp{-t(aA - ei) • (5(c ) - e)} . 

i=[t(aA— ei )] 

A routine application of Borel-Cantelli lemma for integer values of t implies 
the result. □ 

5. The Escape model. 

Proof of Theorem^ (sketch). Since A(t) can not grow faster than the in- 
fected set in the Richardson model with rate 1, from the results of Section 
(use (115(1 ) it follows that for any initial configuration (A(0),B(0)) and all 
sufficiently large n 

E M n (n/c) < exp {— n ■ g(c) + ne} . 

Thus, using the same lines of argument as in Proposition it may be shown 
that for all large t 

A(t) C D c t(ri _ e) , 
A(t) C D t{r2+t) . 
Furthermore, for every c G (?"ij r 2) and any e > 

limsup — log (M n (n/c)) < —g(c) + e. 
n 

(Similarly to the proof of Proposition |2J apply Markov's inequality to get 

P [M n (n/c) > exp {n(-g(c) + e)}] < 
< E [M n (n/c)\ ■ exp {ng(c) — ne} < exp {—ne} 



- log (M n (n/c)) > -g(c)+e 
n 
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and the claim follows by Borel-Cantelli lemma.) Next, we show that on the 
event of survival 



Note that for every non-trivial configuration (A(0), B(0)), for all large t and 
all x G A(t) we have T+(x) D B(t) = 0. Furthermore, almost surely on the 
event of type 1 survival, for every integer m there exist < t < oo and a 
vertex x with m-predecessor y such that x G A(t) and B(t) D T + (y) = 0. 
Define Z\{x) to be a subset of Cij.| +TO n 1+(x) such that ^ G Z\{x) if and 
only if there is a iype 1 poi/« from to + m/c) and there are no 

directed paths from (y,t) to (x,t + m/c). In general, for k > 2, Z/~(x) is 
defined as a subset of vertices in Cui +m ^ n such that G Z^(x) if 

and only if 

(1) 2fc_i G Zk-i(x) where Zk-\ is the m-predecessor of zt; 

(2) there is a iype 1 path from t + (k — l)m/c) to (z^, t + km/c); 

(3) there are no directed paths from (zk-2,t + {k — l)m/c) to (z^-i, t + 
km/c) where Zfc-2 is the m-predecessor of Zk-i ■ 

Let Zk{x) be the cardinality of the set Z^ix). It is clear from the definition 
that (Zk(x))k>o is a Galton- Watson process with the mean offspring number 



Thus, for all sufficiently large m, E[Zi(a;)] > 1. Note that M| x | +m fc(i+fcm/c) 
dominates Zk(x), and hence, on the event of nonextinction of (Zk(x))k>o, 
holds. Observe that for every m, almost surely on the event of survival 
of the type 1, there are infinitely many vertices x at which the Galton- 
Watson processes (Zk{x))k>o can be originated. Hence, (jlfij) holds almost 
surely on the event of nonextinction of the first type. This finishes the proof 



Proof of Theorem^ (sketch). Fix A G (1, A c ), and consider any non-trivial 
initial configuration. With positive probability, there exists a vertex x and a 
Galton- Watson process (Zk(x))k>o (constructed in the proof of Theorem^ 
with EZi(x) > 1. Hence, the Galton- Watson process survives with positive 
probability and so does type 1. 

Consider the case A G (A c , oo). Since A(t) is dominated by the Richardson 
model with rate 1, by Proposition |31 in Section type 1 dies out almost 
surely. □ 



(16) 



liminf — log (M n (n/c)) > —g(c) — e. 



E[Zi(x)] = exp{— mg(c) + o(m)}. 



of Theorem El 



□ 
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